Dynamical properties of the spin-Peierls compound a'— NaV 2 Os 
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Dynamical properties of the novel inorganic spin-Peierls compound a'-NaV20s are investigated 
using a one-dimensional dimerized Heisenberg model. By exact diagonalizations of chains with up 
to 28 sites, supplemented by a finite-size scaling analysis, the dimerization parameter S is determined 
by requiring that the model reproduces the experimentally observed spin gap A. The dynamical 
and static spin structure factors are calculated. As for CuGeOs, the existence of a low energy 
magnon branch separated from the continuum is predicted. The present calculations also suggest 
that a large magnetic Raman scattering intensity should appear above an energy threshold of 1.9 A. 
The predicted photoemission spectrum is qualitatively similar to results for an undimerized chain 
due to the presence of sizable short-range antiferromagnetic correlations. 

PACS numbers: 64.70.Kb, 71.27. +a, 75.10.Jm, 75.40.Mg, 75.50.Ee 
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Recently, the quasi-one-dimensional (ID) compound 
a'-NaV2C>5 has received considerable attention since it 
appears to be the second inorganic material showing a 
spin-Peierls (SP) phase — the first one being CuGeC>3 
jij. Below a transition temperature Tgp « 34AT, the 
compound undergoes a lattice distortion with the open- 
ing of a spin gap. The structure of a'-NaV205 is made 
of quasi-two dimensional layers of VO5 square pyramids 
separated by Na-ions @] . Two types of VO5 chains alter- 
nate : V 4+ Og _ and V^+Og - (V 4+ carries a spin \ while 
V 5+ does not). NaV2C>5 is a good candidate for a ID 
magnetic system since the magnetic V 4+ Og - chains are 
isolated by non magnetic V 5+ Og~ ID structures. 

Originally, the presence of the SP phase transition was 
suggested by experiments on polycrystalline samples [|]]4|] 
which showed a rapid reduction of the magnetic suscep- 
tibility below Tsp — 34 K. From the dependence on the 
orientation of the magnetic field, recent magnetic suscep- 
tibility measurements on single crystals unambiguously 
established the nature of the low-temperature phase 
which is a spin symmetric singlet ground state ||. The 
observation of structural distortions by X-ray diffrac- 
tion ||, NMR Q and Raman scattering || further sug- 
gested that an underlying spin-phonon coupling is re- 
sponsible for the SP transition. Note also that the critical 
temperature T$p is the highest of all known organic or 
inorganic SP compounds (e.g., T$p (CuGeOs) ~ 14 K). 

The magnetic susceptibility x(T) in the high- 
temperature phase above T$p of this compound seems 
well-described by a ID antiferromagnetic (AF) Heisen- 
berg model HQ. Indeed, the crystallographic structure 
of this material suggests that the magnetic frustration is 
very small. Fits of x{T) below T$p provide estimates of 
both the nearest-neighbor spin exchange J and of the 



spin gap A. Recent measurements on single crystals || 
led to J = 441 K and A ~ 85 K in good agreement with 
previous estimates of J [p],[7| and A (4|^|. Then, in our 
analysis below a ratio A/J = 0.193 is assumed. 

Following the approach used for the SP compound 
CuGeC>3, we will consider here a spin-1/2 AF Heisen- 
berg model with an explicit dimerization of the exchange 
coupling to account for the lattice distortion, 



H= J^(l + <*(-!)') 



i+l 



(1) 



Note that in the case of CuGe03 an additional frustra- 
tion was needed to describe the compound (1,^. The 
interchain couplings, although crucial to obtain a finite 
ordering temperature, are expected to be small and will 
be neglected here. 

The model Eq.(l) has a non-zero spin gap for all 8 > 0, 
and first we will determine the value of 8 that reproduces 
the experimentally measured spin gap. The extrapola- 
tion to an infinite chain L — > 00 is performed accurately 
using the scaling law A(L) = A + £ exp(--jt) 0. The 
presence of a spin gap induces a length scale Lq and fi- 
nite size effects are negligible when L ^> Lq. As observed 
in Fig. 0(a) for 8 = 0.05, this scaling behavior is indeed 
accurately satisfied. In addition, for such parameters we 
have found Lq w 18 lattice spacing and, thus, extrap- 
olations using data for systems with up to 28 sites are 
expected to have small error bars. The behavior of the 
spin gap as a function of 8 is shown in Fig. |l|(b). A com- 
parison with the experimental value A/J = 0.193 gives 
an estimate 8 ~ 0.048 for the actual NaV205 compound 
to be compared with 8 ~ 0.014 obtained for CuGe03 ||. 
It is interesting to notice that the dimerization is larger 
for a'-NaV205, although the ratios A/J are similar in 
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both systems pl]| . The reason is that, in contrast to 
a'-NaA^Os, a large frustration exists in CuGeC>3 from a 
sizable next-nearest neighbor coupling constant J'. The 
frustration J' alone can produce a gap when the ratio 
a = J' I "J is larger than a c ~ 0.2411. For CuGe0 3 
a ~ 0.36 was proposed 0, while here a = is assumed. 
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dynamics. Based on our model Eq.(l) and the pa- 
rameter S calculated here, we can predict the dynam- 
ical spin structure factor S zz (q,u>) measured by INS. 
S zz (q,uj) is given by Eq.(||) with A = S z (q) and S z (q) — 
l/^/L^2-exp(iqrj)S z (j). The results on a 28 site chain 
are shown in Fig. ^| for all momenta q — njj, n — 
0, • • • , 14. We clearly observe a well-defined g-dependent 
low energy feature of bandwidth ~ 1.6 J having the 
largest weight located around q = tt. This is certainly 
reminiscent of the Des Cloiseaux-Pearson jl3| (DP) ex- 
citation spectrum of the Heisenberg chain. However, im- 
portant differences arise from the presence of a spin gap: 

(i) there is no intensity for uj < A at q — and q = tt; 

(ii) the lowest singlet-triplet excitation branch which has 
been interpreted as a spinon-spinon bound state fli~4[ is 
well separated from the continuum by a second gap. 
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FIG. 1. (a) Spin gap A(L) in units of J for 8 — 0.05, as a 
function of 1/L. The dashed line is the fitting curve described 
in the text, (b) Bulk extrapolated value of A/ J vs 8. 

Let us now proceed with the study of some dynamical 
properties of a'— NaViOs. Typically, there are a num- 
ber of experiments giving access to frequency dependent 
spectral functions of the form 



I a (uj) = lim 3m (* |^4 — 

7T e^O UJ - 
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where is the ground state, E its energy, and A is 
some operator describing the physical process under con- 
sideration. Using exact diagonalization (ED) techniques, 
Ia(u) can be calculated with a continued fraction expan- 
[12[. An imaginary component ie is added to u> in 
providing a small width to the (5-functions. 
In particular, inelastic neutron scattering (INS) is 
an accurate momentum dependent probe of the spin 




FIG. 2. Spectral function S zz (q,uj) calculated at momenta 
q = ^jp- (L = 28 and 8 = 0.048) and using a small broadening 
e = 0.04J. From bottom to top, n moves from to 14. 

The dispersion relation of the lowest energy magnon 
branch is presented in Fig. with an infinite size extrap- 
olation for momenta q — tt/2 and q = tt (spin gap). The 
finite size effects are quite small, especially for q = tt/2. 
The second peak, as well as the upper limit of the con- 
tinuum of excitations, are also shown. The dispersions of 
the lowest excitations are symmetric with respect to tt/2, 
reflecting the doubling of the unit cell by dimerization. 
However, note that the spectral weight is not symmetric. 
As observed in Fig. ||, we have also explicitly checked at 
q = tt/2 that the magnon excitation is separated from 
the continuum by a gap. 

The static structure factor S zz (q) — J duj S zz (q, uj) and 
the weight of the first peak are shown in Fig. S zz (q) 
is sharply peaked at q = tt due to strong short-range AF 
correlations. A sizable fraction of the weight is located 
above the magnon branch, specially at intermediate mo- 
menta such as q ~ 5tt/7 where the continuum should be 
better observed experimentally. 

It is interesting to make a quantitative comparison 
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FIG. 3. Momentum dependence of the first (o), second 
peak (□) and the upper limit for the continuum excitations 
(o) [L — 28 and 5 = 0.048). The extrapolations to infinite 
size are also shown for q = ir/2 (■, •) and q = n (•). Units 
on the right are meV (assuming J — 440 K). 

with CuGe0 3 . Assuming J = M0K, Fig. || shows 
that the maximum of the magnon branch occurs around 
LOmax — 60 meV while the spin gap is of order A ~ 7.3 
meV. Thus, the energy scales are approximately 4 times 
larger than for CuGeC>3 [^5|Jl^ 1 which might restrict the 
INS experiments to the bottom of the spectrum around 
1 = 7T @- Also, we observed that the ratio uj max /J is 
close to the DP value of 1.57, while in CuGeOa it is ap- 
proximately 1.2 jl^Jl^]. This is due to the fact that the 
large frustration J' in CuGeC>3 affects the entire excita- 
tion spectrum, while in a'-NaV20s only the low energy 
part of the spectrum is modified by the small dimeriza- 
tion. The ratio t-u max /J is then a key quantity to con- 
firm experimentally the absence of frustration in this sys- 
tem using INS. In addition to the change in the value of 
Umax/ J, frustration would also lead to a qualitatively dif- 
ferent global structure of the spectrum. As example, the 
upper limit of the continuum should be better defined for 
CuGe0 3 (J' # 0) than for a'-NaV 2 5 (</' = 0). 

Raman scattering is another powerful technique to 
probe the spin dynamics. The effective Hamiltonian for 
the photon-spin interaction is given by |is| | 



H e s = g 2^ {^in-Rij)(e ou t.Rij)Si.Sj 



(3) 



&in ifiout) is the polarisation vector of the incoming (out- 
going) photons, the sum is over nearest-neighbor spins, 
Rij is the vector connecting them, and g is a coupling 
constant that depends on the incoming photon frequency, 
.ffeff is a spin-singlet, translationally invariant operator, 
and it corresponds to physical processes involving the 
simultaneous excitations of two magnons with opposite 
momenta. Since the small interchain coupling has been 
neglected, Rij has to be collinear to the chain. The 
largest Raman scattering intensity is thus expected for 



FIG. 4. Static structure factor S zz (q) (o) and weight of the 
first excitation (□) as a function of q (L = 28 and S = 0.048). 



a polarisation of both photons along the chain direction. 
The Raman operator can then be written Jl(],[liJ (taking 
g = 1) as Hr = $^(— l)*5i • Si+i . The Raman intensity 
Ir(uj) (Fig. ||) reveals a large scattering band centered at 
a mean energy (defined as the first moment of the spec- 
trum (uj ) = J duLoI R (uj)/ J dwI R (u)) of - 2.9J. I R {u>) 
is fairly smooth (the oscillations at low energy are fi- 
nite size effects) and no Van Hove singularity is observed 
at the energy 2 uj max , associated with the top of the 
magnon branch. Due to the spin gap, an energy thresh- 
old A' appears in the Raman scattering spectrum. The 
infinite size extrapolation of the corresponding singlet- 
singlet gap gives A' ~ 0.37J as indicated in Fig. |. The 
ratio of the singlet-singlet gap over the singlet-triplet gap 
is equal to A'/A ~ 1.9 close to the prediction of V3 Q. 
Since we probe here double magnon excitations this sug- 
gests that magnons are almost non-interacting bosonic 
excitations (which would give A'/A = 2 exactly). 

We end our study of dynamical properties of a'— 
NaV205 with an investigation of the angle resolved pho- 
toemission spectrum (ARPES). In this case the relevant 
operator A in Eq.(^|) is the destruction operator c pa of 
an electron with momentum p and spin a. Results will 
be presented for two hole hopping amplitudes, i.e. t = J 
and t = 2.5 J, since it is a priori difficult to anticipate its 
actual value. The results shown in Fig. ^ are very sim- 
ilar to the case of a single hole in a half-filled infinite - 
U Hubbard model |20| or in the frustrated Heisenberg 
chain The overall scale of the spectrum is clearly 

given by t. For p < tt/2 a holon and spinon branches 
appear, while a "shadow band" is observed for p > tt/ 2 
caused by short-range magnetic scattering at q — tt [[15| . 
Above the continuum, the high energy structure can be 
associated with a reflection of the shadow band at the 
zone boundary p0| |. However, note that the presence of 
a spin gap introduces some subtle differences with pre- 
dictions for undimerized systems that could be detected 
in ARPES experiments. For instance, doping a spin gap 
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FIG. 5. Raman intensity for L = 28 and 6 = 0.048. A 

broadening e = 0.2 J was used. The arrows indicate the ex- FIQ 6 Hole gpectral function A ^ w) on a chain of 20 sites 

trapolated singlet-singlet gap A at low energy and the first for t/J = 1 (a) afld f/J = 2 g and broadening e = .1J. 

moment of the distribution (w) at higher energy. FrQm bottQm tQ tQp the momentum varies from to 7r . 



insulator leads in general to a metallic state with only one 
zero-energy mode ||l[] (corresponding to collective charge 
excitations only) . This is reflected in Fig. || by the fact 
that the so-called spinon branch is in fact a broad struc- 
ture instead of a branch-cut. In addition, the small peak 
at very low energy just above kp — n/2 might be associ- 
ated to a holon-spinon bound state. Similar results were 
obtained for CuGe0 3 ]l5p . In spite of these subtleties, it 
is clear that the spectral function of the dimerized model 
Eq.(l) has strong similarities with undimerized systems. 
The main reason is that here an explicit dimerization co- 
exists with sizable short-range AF correlations, a detail 
not sufficiently remarked in the literature on the subject. 

In conclusion, using recent experimental data on the 
SP a'-NaV2 05 system, the magnitude of the dimeriza- 
tion of the AF exchange coupling along the chain has 
been determined. In the framework of a ID dimerized 
Heisenberg model, several theoretical predictions for INS, 
Raman double magnon scattering, and ARPES were here 
presented. Our calculations are expected to provide the- 
oretical guidance for future experiments on SP systems. 
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